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Abstract  --  Spatially  resolved  electron  temperatures,  electron  number 

densities  and  gas -kinetic  temperatures  have  been  determined  at  five 

observation  heights,  nine  radial  positions,  and  two  settings  of  RF  input 

power  in  an  aerosol- free  argon  ICP  sustained  with  a  commercial  MAK  torch. 

measurements  were  obtained  by  means  of  simultaneous  ruby- laser  Thomson 

scattering  and  Rayleigh  scattering.  In  the  central  channel  of  the  ICP  the 

electron  temperature  ranges  from  7000  K  to  8500  K,  the  gas  kinetic 

temperature  is  about  1000  K  to  2000  K  lower  than  the  electron  temperature, 

15  -  3 

and  the  electron  number  density  varies  from  1.9  to  7.0  x  10  cm  .  Among 
the  measured  parameters,  the  electron  number  density  is  most  sensitive  to 
the  input  RF  power  and  is  noticeably  higher  than  the  LTE  value  determined 
from  the  measured  electron  temperature  and  the  Saha  equation.  The 
overpopulation  of  electrons  can  probably  be  attributed  to  fast  ambipolar 
diffusion  and  slow  ion-electron  recombination  processes,  resulting  in  a 
recombining  plasma  in  the  analytical  zone.  This  non-LTE  feature  might  be 
significant  for  excitation  and  ionization  mechanisms  in  the  ICP. 
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1.  INTRODUCTION 

If  the  inductively  coupLed  argon  plasma  (ICAP)  sustained  at  atmospheric 
pressure  were  in  thermodynamic  equilibrium  (TE) ,  a  single  physical  parameter 
(e.g.  temperature  T  or  electron  number  density  ng)  would  be  sufficient  to 
completely  describe  its  behavior.  That  is,  under  TE  conditions,  electron 
number  density,  the  number  density  of  ions  at  different  ionization  stages, 
energy- level  populations,  the  frequency  distribution  of  radiation  and  the 
velocity  distribution  of  all  species  could  be  all  exactly  determined  from 
the  ideal  gas  law,  the  Boltzmann  distribution,  the  Saha  equation,  Planck's 
radiation  law  and  the  Maxwell  distribution. 

If  the  ICAP  is  not  in  TE,  but  were  instead  in  local  thermodynamic 
equilibrium  (LTE) ,  it  could  again  be  described  by  only  one  physical 
parameter.  However,  the  parameter  would  be  a  function  of  spatial  position 
e.g.  T(r,z)  or  ne(r,z),  where  r  and  z  are  radial  position  and  observation 
height,  respectively.  For  a  given  spatial  point  (r,z),  all  the  detailed 
information  mentioned  above  (except  the  Planck  radiation  law)  could  be 
determined  by  T(r,z)  or  ne(r,z),  using  the  same  equations  or  laws  as  for  the 
case  of  TE.  Unfortunately,  as  pointed  out  quite  early  by  many  investigators 
[1-5],  the  ICAP  is  not  in  LTE  either.  Thus,  even  at  a  given  spatial  point, 
a  single  physical  parameter  is  not  sufficient  to  describe  the  plasma 
bena/ier.  Of  course,  analyte  excitation  might  be  attributable  to  the 
presence  or  concentration  of  particular  species  (e.g.  Ar  or  Ar+) .  However, 
measuring  the  concentrations  of  such  species  alone  does  not  characterize  the 
plasma  as  a  whole  and  might  not  e”on  elucidate  ^n«ly to -excitation  events. 

The  lev  quest  ion  then  is.  "What,  are  the  most  useful  fundamental  parameters 
:  r , :  .I,!  ;  in--  th.-  b.  ii  i  v ;  or  :  n  ■  I' •A;’  "•  .  Tin  nnswr  i  ■.  loe.il 
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The  reasons  for  this  conclusion  are  compelling.  When  the  sample 

aerosol  is  introduced  into  a  plasma  through  the  central  channel,  the  first 

events  involved  are  desolvation,  vaporization  and  dissociation,  which  are 

greatly  affected  by  the  gas  temperature  T  .  Electron  temperature  T  and 

g  e 

electron  number  density  ng ,  on  the  other  hand,  are  the  determining 

parameters  for  electrical  conductivity  and  ambipolar  diffusion  in  the 

plasma.  Two  other  important  transport  coefficients  -  thermal  conductivity 

and  viscosity  -  are  controlled  by  T  n  and  T  [6]. 

e  e  g 

Furthermore,  Tg  and  ng  are  extremely  important  in  evaluating  a  variety 
of  plausible  analyte  excitation  and  ionization  mechanisms  in  the  ICAP. 
First,  the  rate  constants  for  collisional  excitation,  de-excitation, 
ionization  and  three-body  recombination  are  generally  sensitive  to  Tg ,  and 
are  all  proportional  to  ng  (except  for  three -body  recombination  which  is 
proportional  to  the  square  of  ng  [7-9]).  There  is  little  doubt  that  when 
both  Tg  and  ng  are  high,  electron  impact  will  dictate  the  energy- level 
populations  of  analyte  atoms  and  ions,  and  should  play  a  vital  part  in  the 
excitation  and  ionization  processes.  Of  course,  the  Boltzmann  distribution 
and  the  Saha  equation  may  not  apply  to  such  events,  as  ng  might  not  be 
the  LTE  value  corresponding  to  Tg .  Secondly,  Penning  ionization  would 
contribute  to  analyte  ionization  and  excitation  only  when  argon  metastable 
levels  are  highly  populated  [2];  this  event  would  occur  only  when  electrons 
are  overpopulated  with  respect  to  Te  [10],  Thirdly,  charge  transfer,  which 
Has  rec^nrlv  been  suggested  to  play  a  significant  role  in  the  ICP  [11,12 j. 
is  indirectly  influenced  by  T  and  n  If  T  is  quite  high,  manv  element'- 
••••<■>';!  d  be  ionized  hv  .-•!  r>,j  :  i:- - 1  ,  and  m.lv  a  sn  all  port  ion  •>:  ' 


ana  1  y  to  spec  i  os  would  remn  i  n  as  atom:;.  Thu:; ,  chart’,*’  t  ra:i:;l  .*r  would  not 
contribute  apprec  i  ah  Lv  to  the  degree  of  ioniznt  ion.  On  the  other  hand,  if 
electrons  are  overpopulated  and  T  remains  relatively  low,  charge  transfer 
might  he  of  importance. 

From  the  above  discussion,  local  T  n  and  T  values  are  critical  in 

e  e  g 

describing  the  macroscopic  behavior  and  microscopic  processes  that  occur  in 
the  ICP.  It  is  therefore  not  surprising  that  many  authors  have  measured 
these  parameters  using  different  methods  of  plasma  diagnostics  [1,3-5,13- 
35],  Among  these  methods,  laser  light  scattering  offers  a  number  of  unique 
advantages  [36-39],  Most  importantly,  when  a  focused  laser  beam  is 
scattered  by  the  plasma,  both  Thomson  scattering  and  Rayleigh  scattering  can 
be  measured  at  the  same  time,  making  possible  the  simultaneous  measurement 
of  T  n  and  T  on  a  spatially  and  temporally  resolved  basis.  These 

e  e  g 

procedures  greatly  reduce  the  errors  generated  by  the  Abel  inversion  required 
by  most  other  methods  and  avoid  the  plasma  drift  (instability)  problem,  at 
least  for  any  chosen  spatial  point  of  the  measurements. 

In  the  present  paper,  we  report  local  values  of  T  ,  n  and  T  measured 

e  e  g 

simultaneously  by  means  of  Thomson  scattering  and  Rayleigh  scattering  from 
an  aerosol-freelCAP  sustained  with  a  commercial  MAK  torch.  The  measurements 
were  performed  at  five  observation  heights  (5.0,  7.5,  10.0,  12.5,  and  15.0 
mm  above  the  load  coil  (ALC))  and  at  nine  radial  positions  (-6.0,  -4.5,  - 
3.0,  -1.5,  0,  1.5,  3.0,  4.5  and  6.0  mm  from  the  plasma  axis).  Two  settings 
of  RF  input  power  (0.75  and  1.0  kW)  were  used  for  comparison.  Compared  to 
the  results  of  a  previous  study  [28],  the  Tg  and  ng  values  presented  here 
should  be  more  reliable,  because  a  new  data  treatment  procedure  has  been 
section).  In  addit  ion. 
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i.K’iil  pus  lav  and  Rayleigh  scattc  r  ini’  has  Iici-ii  m  f  n  i  n;  i  :-:<*<  1 .  Thi'nivt  ir;i  I 
formulas  employed  in  the  calculat  ions  of  T  ,  n  and  T  ,  f  rom  the  scattered 
light  signals  will  be  briefly  reviewed,  and  some  new  calibration  and 
normalization  procedures  for  the  multichannel  Thomson  scattering  system  will 
be  described.  Finally,  the  experimental  results  will  be  discussed  and 
conclusions  will  be  drawn. 
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2  .  THEORY 

In  this  snot  ion,  the  theoretical  formulas  used  for  c;i  1  on  i  a  t  i  m;  T\  ,  n 
and  T  from  tin-  measured  Thomson-  scat  tor  i  ng  and  Rayleigh-scat  t  oring  signals 
will  be  briefly  summarized.  For  a  more  detailed  discussion  of  the  theory 
from  which  the  formulas  are  derived,  the  reader  is  referred  to  earlier 
papers  [35,39] . 

A  Thomson- scatter ing  signal  1(A),  expressed  as  a  function  of  wavelength 
A  is  related  to  the  shift  AA  from  the  incident  laser  wavelength  AQ  by  the 
formula 

me  °2  2 

In  1(A)  =  C - - - - -  ‘  (AA)  ^  (1) 

8kTesinz(0/2)‘  a£ 

where  C  is  a  constant  that  depends  on  the  entire  scattering  system,  mg  is 

the  mass  of  an  electron,  c  the  speed  of  light,  k  the  Boltzmann  constant,  and 

6  the  angle  between  the  observation  direction  and  the  laser  incident  beam. 

According  to  Eqn  (1) ,  the  local  electron  temperature  Tg  can  be  determined 

from  the  slope  of  a  plot  of  lnl(A)  vs.  (AA)2.  To  use  Eqn  (1),  one  must  be 
2 

sure  that  mec  and  kTg  are  in  the  same  units,  no  matter  whether  Tg  is  in  eV 

O  p 

or  in  K.  The  same  is  true  of  course  for  (AA)  and  A  Q. 

Because  the  amplitude  of  the  entire  Thomson- scattering  spectrum  is 
proportional  to  the  electron  number  density,  the  sum  of  the  absolutely 
calibrated  signals  from  all  spectral  channels  can  be  used  to  calculate  the 
local  value  of  ng  as  shown  in  Eqn  (2) . 

J  1(A) dA 


"ra 

( 2 ) 
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'.ihrro  n.  is  the  number  density  of  argon  atoms  tit  the  tempo  rat  ure  tit.  which 
the  Rayleigh-scattering  measurement  is  performed  (usutilly  room  temperature 
[35,39]),  the  integration  term  expresses  the  sum  of  Thomson- scatter ing 
signals  from  all  spectral  channels,  1^  is  the  Rayle igh - scattering  signal 
from  the  central  spectral  channel  located  at  the  incident  laser  wavelength, 
is  the  Rayleigh-scattering  cross  section  of  argon  gas,  and  a T  is  the 
Thomson- scattering  cross  section.  In  Eqn  (2),  Rayleigh  scattering  serves  to 
calibrate  on  an  absolute  basis  the  Thomson  scattering  signals  [35,39], 

When  the  plasma  is  on,  the  Rayleigh-scattering  signal  will  be  greatly 
reduced  compared  to  that  when  the  plasma  is  off.  This  is  because  the 
magnitude  of  the  Rayleigh  signal  is  proportional  to  the  argon-atom  number 
density  which,  in  turn,  is  inversely  proportional  to  the  gas  temperature 
according  to  the  ideal  gas  law.  Therefore,  the  ratio  of  the  Rayleigh- 
scattering  signal  when  the  plasma  is  off  IpA(off)  to  that  when  the  plasma  is 
on  IpA(on)  can  use<*  to  determine  the  local  gas  temperature  according  to 
Eqn  (3) . 


IRA(of f ) 

T  =  T 

g  *  room 

^RA(on) 


(3) 


where  Troom  is  the  room  temperature  at  which  is  measured.  Note 

that  in  Eqn  (3),  both  IpA(off)  and  IpA(on)  should  be  pure  Rayleigh- 
scattering  signals;  in  other  words,  the  contribution  of  stray  light  to  the 
central  -  spec tral -channe 1  signal  must  be  subtracted.  It  should  be  pointed 
out  also  that  Eqns.  (1)  and  (2)  for  the  determination  of  T  and  n  are  valid 
ou!v  when  i  hi’  so -ra  1  ]  ■■•d  a  !  f  or  i  u  >•. ,  lam- • ;  .  ■  r  n  '  s  much  smaller  than  ’  ' 


1 


.>'*  '■  .  <»:  horvi  s<*  a  look-up  t  ah  1  e  can  he  used  to  1  i  nd  the  true  T  and  n  1  rum 

t  !.>•  values  of  T  and  n  calculated  bv  usinr  t  hose  two  pquation;:.  For  m* » l'-  • 
detail  about  this  data- treatment  procedure,  the  reader  is  referred  to  a 
recent  paper  [40]. 


3 .  EXPERIMENTAL 

The  entire  multichannel  Thomson- scattering  system  constructed  in  our 
laboratory  has  been  described  previously  [41].  In  the  present  study,  a 
commercial  MAK  torch  with  relatively  low  argon  gas  flow  and  RF  input  power 
was  used.  Two  settings  of  ICP  operating  conditions  used  in  the  measurements 
are  listed  in  Table  1.  For  these  conditions,  the  tip  of  the  "bullet"  in  the 
ICP  is  located  at  2  mm  ALC  (0.75  kW)  and  0  mm  ALC  (1.0  kW) .  In  order  to 
improve  accuracy,  the  multichannel  detector  system  was  normalized  by  using  a 
calibrated  tungsten  ribbon  lamp,  located  exactly  in  place  of  the  plasma.  In 
earlier  work  [28,41],  a  light-emitting  diode  (LED)  was  mounted  on  the  top  of 
the  focusing  concave  mirror  inside  the  monochromator.  Light  emitted  by  the 
LED  was  used  to  illuminate  directly  the  whole  area  of  the  fiber-optic  array 
but  without  passing  through  any  of  the  optical  elements  between  the  ICAP  and 
the  monochromator.  In  contrast,  the  light  emitted  by  the  tungsten  ribbon 
lamp  travels  exactly  the  same  path  and  illuminates  the  same  portion  of  the 
fiber-optic  array  as  does  the  scattered  laser  light.  In  addition,  the 
illumination  by  the  tungsten  ribbon  lamp  also  simulates  a  scattering  signal 
well  because  it  is  pulsed  by  rotation  of  a  mirror  located  in  the  light¬ 
collecting  system.  Of  course  the  pulse  length  (about  25  ns)  is  much  longer 
than  that  from  the  earlier  LED  (400  ns).  During  the  experiment,  the  LED  was 
still  used  to  monitor  the  stability  oi  each  spectral  channel.  To  compare 
two  approaches.  t.ho  norma  1  i /.a  *  ion  factors 
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for  the  first  nine  spectral 
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Brc.-iuso  the  width  ami  height  oj  the  monorh  roma  t  o  r  out  ranee  .slit  were 

.set  til  0.8  and  1.0  inm ,  respectively,  and  the  optical  collecting  system  has 

an  image  magnification  of  1:1,  the  spatial  resolution  in  our  measurement  was 
3 

better  than  1  mm  .  The  temporal  resolution  is  limited  by  the  pulse  length 
of  the  Q-switched  ruby  laser  (-  25  ns). 

The  stray- light  contribution  to  the  signal  from  the  central  spectral 
channel  was  determined  by  measuring  the  Rayleigh- scatter ing  signals  from 
argon  and  helium,  respectively;  both  gases  were  at  atmospheric  pressure  and 
at  room  temperature  [35], 

4 .  RESULTS 

Radial  profiles  of  T  n  and  T  at  five  observation  heights  (5.0,  7.5, 

e  e  s 

10.0,  12.5  and  15.0  mm  ALC)  are  shown  in  Figures  1-3  for  0.75  kW  of  RF  input 
power.  Similar  profiles  for  1.0  kW  of  RF  input  power  are  presented  in 
Figures  4-6.  In  general,  at  each  observation  height,  nine  radial  positions 
including  the  spatial  point  on  the  plasma  axis  were  chosen  with  an  increment 
between  adjacent  locations  of  1.5  mm.  Unfortunately,  the  electron  number 
density  drops  rapidly  wTith  increasing  observation  height,  particularly  at 
larger  radial  distances  from  the  plasma  axis  (see  later  discussion), 
resulting  in  very  noisy  Thomson  scattering  signals  at  the  plasma  edge. 

These  noisy  signals  arise  not  only  because  the  Thomson  scattering  signal  is 
proportional  to  n  ,  but  also  because  stray  light  remains  almost  unchanged  at. 
the  plasma  edge,  and  the  Rayleigh-scattering  signal  becomes  even  higher. 

This  latter  increase  is  flue  to  the  cooler  gas  temperature  at  the  plasma 
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■  :  I:1.!)  m:n  thoro  are  j  ei.it  .  >  s!iuv;n  in  I-'i.i;urcs  4-6  (1  .0  im: 

.  .  i  <  .  i  t  .  t  in  Figues  1-  i  (()./')  ill.')  ,  because  the  elect  ron  1 1  n ; : . 1  >  *  •  l  ('.•■unity 
:  teases  rapidly  with  RF  input  power  (sue  later  discussion)  . 

It.  is  appropriate  to  question  whether  the  noisty  Thomson  signals  at  the 
plasma  edge  can  be  improved.  The  Thomson-scattering  signals  themselves  are 
currently  photon-noise  limited.  As  a  result,  they  can  be  improved  only  be 
acquiring  larger  numbers  of  photons.  Importantly,  this  increase  is  not 
possible  through  use  of  a  higher-power  source  laser,  since  a  higher 
irradiance  would  perturb  the  plasma.  Instead,  it  would  be  desirable  to 
integrate  a  larger  number  of  laser  shots,  possilby  either  by  increasing  the 
signal- integration  interval  or  by  using  a  laser  with  a  higher  pulse 
repetition  rate.  The  latter  strategy  is,  of  course,  the  most  attractive. 

The  background  noise  at  the  plasma  edge  is  greater  than  that  in  the  plasma 
center,  a  fact  which  further  degrades  Thomson  signals  acquired  there. 
Fortunately,  the  contribution  of  this  noise  to  a  multichannel  Thomson 
spectrum  arises  as  a  shot-noise-limited  background  and,  like  the  Thomson 
signal  itself,  can  be  improved  by  accumulating  more  laser  shots. 

From  Figures  1  and  4,  the  radial  distributions  of  Te  are  fairly  flat 
below  10  mm  ALC  for  both  RF  input  power  settings but  become  more  bell- 
shaped  at  higher  observation  positions  (cf.  Figures  1  and  4).  In  contrast, 

T  (Figures  3  and  6)  decreases  rapidly  with  increasing  radial  distance  from 

o 

the  plasma  axis  even  at  low  observation  heights.  Furthermore,  T  is  usually 

O 

about  1000  K  to  2000  K  lower  than  Te  at  every  plasma  location.  It  is  also 
apparent  that  the  radial  profiles  of  ne  are  quite  symmetric  with  respect  to 
the  plasma  axis  at  almost  all  observation  positions  (cf.  Figures  2  and  5). 
The  value  of  n  drops  more  rapidly  with  increasing  radial  distance  from  the 
:  mi'  i  .  i i  when  :  he  oh.se  rvat  ion  te-  i  i'.ht  becomes  larger.  None  of  the  radial 
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!  i:  ;  If  ch.ii'.i’.i'  in  1'  is  ‘<>\uui  at.  olisci'v.it  ion  heights  lower  than  10  n>:n  ALC. 
Heva-ver,  an  apparent  iunvasc  (about  1)00  K)  in  T  occurs  at  higher 
observation  positions.  By  comparison,  the  effect  of  RF  input  power  on  T 

o 

can  be  clearly  seen  even  at  5  mm  and  7.5  mm  ALC ,  particularly  at  the  plasma 
edge.  An  obvious  increase  of  n£  with  RF  input  power  can  be  seen  at  all 
spatial  points,  an  effect  which  becomes  even  more  significant  higher  in  the 
plasma.  A  more  detailed  and  quantitative  discussion  of  these  observations 
is  given  in  the  following  section. 


5.  DISCUSSION  AND  CONCLUSIONS 

Electron  temperatures  Tg  as  a  function  of  observation  height  on  the 
plasma  axis  are  shown  in  Figure  7  for  both  settings  of  RF  input  power  (open 
triangles:  0.75  kW,  filled  triangles:  1.0  kW) .  It  can  be  seen  that  Tg 

does  not  change  greatly  with  observation  height  from  7.5  mm  to  15.0  mm  ALC, 
and  that  increasing  RF  input  power  from  0.75  kW  to  1.0  kW  raises  only 
slightly  (less  than  10%)  the  Tg  value  on  the  plasma  axis. 

Similar  plots  for  gas  temperatures  T  are  shown  in  Figure  8  (filled 

O 

triangles:  0.75  kW,  open  triangles:  1.0  kW) .  Compared  to  Te ,  the  Tg  value 

on  the  plasma  axis  decreases  more  quickly  with  observation  height  for  both 
settings  of  RF  input  power.  Also,  the  increment  in  Tg  when  RF  input  power 
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.'u,  s  T  (  Fi  r.urv  /)  or  T  (Figure  S').  (Not*-  {Ik-  1  o;’,ar  i  t  lnni  c  w  ri  i  ca  1  axis  os. 
Figure  11 )  .  1)'  nddit  ion,  n  is  much  more  sons  it  i\v  to  the  RF  input  power 

;  nan  either  'l'o  or  T  ,  particularly  i  n  t  ho  upper  observation  regions.  For 
example,  at  1.0  kW  the  no  value  at  15  mm  ALC  is  about  twice  that  at  0.75  kW . 

Vertical  distributions  of  ng  at  3  mm  off  the  plasma  axis  are  plotted  in 
Figure  10  for  both  settings  of  RF  input  power  (filled  circles:  0.75  kVJ , 
open  circles:  1.0  kW)  .  Obviously,  all  ng  values  at  3  mm  off  the  plasma 
axis  are  lower  than  their  corresponding  ones  on  the  plasma  axis  (Figure  9) . 
Furthermore,  at  3  mm  off  the  plasma  axis,  ng  declines  with  observation 
height  more  rapidly  than  occurs  on  the  plasma  axis  (compare  Figures  9  and 
10),  particularly  at  lower  RF  input  power. 

Now  let  us  compare  in  more  detail  the  electron  temperatures  and  their 
corresponding  gas  temperatures,  and  see  how  different  they  are  from  each 
other.  From  Figures  7  and  8,  the  disparity  between  T  and  Tg  exists  at  all 

C 

observation  heights,  and  is  generally  between  1000  K  and  2000  K,  as 

mentioned  in  section  4.  Ncte  also  that  a  rise  in  RF  input  power  reduces  the 

deviation  slightly. 

A  comparison  of  the  vertical  distribution  of  Te  and  Tg  at  1.0  kW  of  RF 

input  power  and  at  3  mm  off  the  plasma  axis  is  shown  in  Figure  11.  It  is 

not  surprising  that  the  off-axis  difference  between  T  and  Tg  is  greater 

O 

than  that  on-axis  (See  Figures  7  and  8). 

Another  way  to  evaluate  how  far  the  plasma  deviates  from  LTE  is  to 
compare  the  measured  ne  with  values  calculated  according  to  the  Saha 
equation.  Examples  of  such  comparisons  are  given  in  Figures  12  and  13  for 
0.75  kW  and  1.0  kW  of  RF  input  power,  respectively.  In  Figure  12,  the 
filled  circles  are  n  data  measured  on  the  pi  asm a  axis  directly  by  Thomson 
scattering,  whereas  the  open  circle  •  i  r  <  •  values  e  a  I  r  n  1  a  t  <■<!  from  the  Saha 
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apparent  deviation  or  from  the  LTE  valuers  (derived  from  T,)  exists  .at  all 
observation  heights.  In  general,  this  deviation  ranges  from  a  factor  of  3 
to  one  order  of  magnitude,  depending  on  the  spatial  position  and  RF  input 
power.  Thus,  electrons  in  the  plasma  are  overpopulated  compared  to  what  Te 
would  predict,  resulting  in  what  can  be  viewed  as  a  recombining  plasma. 

The  values  of  T  n  and  T  presented  in  Figures  1-6  are  comparable  to 
e  e  g 

those  reported  previously  by  other  authors  [1,3-5,13-35].  Although  the 
electron  temperatures  measured  here  are  about  20%  lower  than  those  reported 
by  Batal  et  al.  [13,14],  their  observation  location  was  lower  (2  mm  ALC)  and 
their  RF  input  power  was  higher  (1.5  kW)  than  ours.  Better  agreement 
(within  10%)  has  been  found  between  our  Te  data  and  those  recently  reported 
by  Davis  and  Du  [16],  All  these  authors  determined  Tg  by  measuring  the 
intensities  of  the  argon  emission  line  and  the  background  at  430  nm.  Lu  et 
al.[17]  more  recently  published  Tg  values  measured  by  Langmuir  probes  in  the 
ICP.  Their  results  are  found  very  close  to  ours  with  the  same  RF  input 
power . 

The  electron  number  densities  shown  in  Figures  2  and  5  range  from  1  x 
10^  to  7  x  1015  cra*^,  which  are  similar  or  comparable  to  those  previously 
reported  [5,13-27],  except  results  from  Langmuir -probe  measurements,  which 
seem  somewhat  too  high  [17].  In  general,  our  ng  values  are  in  better 
agreement  with  those  derived  from  absolute  continuum  intensities  [1,19],  but 
somewhat  higher  than  those  measured  by  line  broadening  at  low  observation 
positions  [25,27],  The  discrepancy  might  be  due  to  the  different  torches 
and  ICP  operating  conditions. 
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T:i«‘  1 ei:  pe  natures  oiled  in  Figures  3  ami  (>  are  in  the  range  reported 

be  u;  In  !  anther::  [  1  ,  1  8 , 3  9 ,  3 1  -  34  j  .  Interest  inglv,  our  data  for  '1'  are  in 
good  agreement  with  those  measured  by  Doppler  broadening  of  emission  lines 
[29,33,34],  but  are  much  higher  than  values  derived  from  molecular 
rotational  spectra  [1,18,31,32].  This  dissimilarity  is  not  surprising, 
since  our  determination  of  T  from  Rayleigh  scattering  is  based  on  the  ideal 

O 

gas  law,  which  relates  the  gas  pressure  with  the  kinetic  temperature  and 
argon-atom  number  density  in  the  plasma.  Because  the  gas  pressure  (and 
therefore  density)  is  produced  by  the  translational  moment  of  the  plasma 
species,  the  T  determined  by  Rayleigh  scattering  should  be  closer  to  the 

O 

Doppler  temperature  than  to  the  rotational  temperature  under  non-LTE 
conditions.  The  lower  rotational  temperature  suggests  also  that  the 
transfer  of  particle  kinetic  energy  to  molecular  rotational  movement  is  a 
slower  process  in  the  ICP  than  is  kinetic-energy  transfer  from  electrons  to 
argon  atoms . 

As  mentioned  above,  the  measured  ng  is  higher  than  the  LTE  value  for  ng 

calculated  from  the  measured  Te .  Viewed  differently,  if  one  uses  the 

measured  ng  to  determine  Tg  according  to  the  Saha  equation  [42] ,  the  Tg  would 

be  overestimated.  This  fact  can  explain  the  recent  findings  of  Blades  et 

al.  [42]  about  the  deviations  of  the  ICP  from  LTE,  when  they  observed  ion- 

atom  emission- intensity  ratios  for  Sr,  Ca,  Mg,  Cd  and  Zn.  They  found  that 

the  experimental  ratios  (I^/Ia)exp  are  several  times  lower  than  the  LTE 

ratios  kW  °f  RF  input  power.  The  reason  is  that  they  used 

their  measurer  nc  to  determine  Te  by  assuming  LTE  conditions  and  then  used 

those  Tg  and  ng  values  to  calculate  (I^/I^^xp-  However,  as  verified  here, 

the  actual  T  is  lower  than  the  T  value  determined  in  their  way,  and  the 
e  c 

ionic  levels  are  actual  lv  1  t  p< pti  1  a  t  ■  d  .  Therefore,  the  measured  ion-atom 
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i  • i  ss  i  0:1  int  ensi  t  y  ratios  should  be  lowiM'  than  tin-  >  i  Tf  1  i  ns’.  "I-T:." 

v.;  i  n.'s  .  Th>‘  same  reason  i  ::g  ]  e.uis  to  their  cone  !  u:  i  on  that  i  tie  ex  |  >e  r  i  mu  tit  a  1 
degree  of  ionization  for  the  same  elements  is  lower  than  that  under  " LT K " 
conditions.  Clearly,  both  Te  and  ng  should  be  measured  independently  but 
directly  in  order  to  better  describe  the  excitation  and  ionization  processes 
in  the  ICP. 

The  overpopulation  of  electrons  in  the  analytical  zone  might  be  a 
significant  non-LTE  feature  for  analytes  in  the  ICP,  because  it  not  only 
favors  the  excitation  and  ionization  by  direct  electron  impact  and  charge 
transfer,  but  also  gives  rise  to  overpopulated  argon  excited  states  and 
high-energy  electrons  through  recombination  and  superelastic  collision  [43]. 
There  are  two  possible  reasons  why  electrons  in  the  ICP  are  overpopulated. 
First,  the  ambipolar-dif fusion  process  [44],  which  causes  electron-ion  pairs  to 
migrate  from  the  fireball  to  the  cooler  plasma  regions,  is  faster  than  the 
thermal -conduction  process  that  leads  to  the  transfer  of  kinetic  energy  from 
the  fireball  to  the  cooler  plasma  regions.  As  a  consequence,  the  analytical 
zone  gains  electrons  faster  than  it  Is  heated  thermally,  resulting  in  an 
electron  number  density  that  is  higher  than  the  LTE  value.  Secondly,  when 
the  plasma  travels  upward,  the  argon  ion-electron  recombination  process  is 
relatively  slow  [45],  compared  to  the  plasma  cooling  caused  by  energy  loss 
through  thermal  conduction,  convection  and  radiation.  Thus,  electrons 
remain  overpopulated  and  the  plasma  remains  in  a  recombining  mode. 

Finally,  the  fact  that  none  of  the  radial  profiles  of  T  n  and  T 
shown  in  Figures  1-6  are  doughnut  -  shaped  is  probably  due  to  the  ICP 
operating  conditions  suggested  by  the  manufacturer  of  the  MAK  torch  (see 
Table  1).  Because  the  central  argon  gas  flow  was  only  0.5  1/min.  and  no 
aerosol  war;  added,  all  the  above -men;  ioned  parameters  would  quickly  peak  ><:' 
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;>1. i  nxi:;.  As  a  matter  of  fact,  our  observation  region  was  in  the 

!  anal  v;  >  ca  1  zone .  whereas  the  initial  radiation  zone  was  found  far 

lower  in  the  plasma  (no  sample  aerosol  was  used  in  this  study). 

We  conclude  that,  under  our  experimental  conditions,  (1)  T  in  the  ICP 

& 

is  about  1000  to  2000  K  lower  than  T  ,  and  ng  is  overpopulated  by  a  factor 

of  3  —  compared  to  the  LTE  value  determined  by  Tg ;  (2)  ne  is  most  sensitive 

of  the  three  variables  to  RF  input  power  and  to  spatial  position;  (3)  fast 

ambipolar  diffusion  and  slow  argon  ion-electron  recombination  probably  lead 

to  the  overpopulation  of  electrons,  resulting  in  a  recombining  plasma  in  the 

analytical  zone;  (4)  the  ambipolar-dif fusion  and  thermal -conduction 

processes  are  fast  enough  to  cause  rapid  plasma  decay;  (5)  T  ,  n  and  T  are 

e  e  g 

the  most  important  parameters  to  describe  the  non-LTE  features  of  the  ICP, 
and  all  these  parameters  can  be  simultaneously  determined  by  Thomson 
scattering  and  Rayleigh  scattering  on  a  spatially  and  temporally  resolved 
basis . 
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ble  L  -  IGF  operating  conditions 


Torch 

RF  Frequency 
RF  input  power 
Outer  argon  flow 
Intermediate  argon  flow 
Central  argon  flow 
Sample  flow  rate 


Commercial 
27.12  MHz 
0.75  kW 
10.0  1/min 
0.5  1/min 
0.5  1/min 
none 


MAK  torch 

1.0  kW 
11.0  1/min 
0.5  1/min 
0.5  1/min 
none 


lable  2.  Comparison  of  Multichannel  Normalization  Factors  Determir 
Tungsten  Ribbon  Lamp  (W)  and  Light  Emitting  Diode  (LED) 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 
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Figure  Captions 

Radial  distributions  o)  electron  temporal  ures  measured  by 
Thomson  scattering  at  five  observation  heights  above  the 
load  coil  and  at  an  RF  input  power  of  0.75  kW. 

Radial  distributions  of  electron  number  densities 
measured  by  Thomson  scattering  at  five  observation 
heights  above  the  load  coil  and  at  an  RF  input  power  of 
0.75  kW. 

Radial  distributions  of  gas  kinetic  temperatures 
measured  by  Rayleigh  scattering  at  five  observation 
heights  above  the  load  coil  and  at  an  RF  input  power  of 
0.75  kW. 

Radial  distributions  of  electron  temperatures  measured 
by  Thomson  scattering  at  five  observation  heights  above 
the  load  coil  and  at  an  RF  input  power  of  1.0  kW. 

Radial  distributions  of  electron  number  densities 
measured  by  Thomson  scattering  at  five  observation 
heights  above  the  load  ceil  and  at  an  RF  input  power  of 
1.0  kW. 

Radial  distributions  of  gas  kinetic  temperatures 


measured  by  Rayleigh  scattering  at  five  observation 


i  i  ■  ■  :  r.ii;  s  ,ibovi‘  t  h<‘  load  coil  and  at  an  RF  input  power  of 
i.a  k\: . 


Figure  /.  Electron  temperatures  on  the  plasma  axis  tit  five 

observation  heights  above  the  load  coil  and  at  an  RF 
input  power  of  0.75  kW  (open  triangles)  and  1.0  kW 
(filled  triangles). 

Figure  8.  Gas  temperatures  on  the  plasma  axis  at  five  observation 
heights  above  the  load  coil  and  at  an  RF  input  power  of 
0.75  kW  (filled  triangles)  and  1.0  kW  (open  triangles). 

Figure  9.  Electron  number  densities  on  the  plasma  axis  at  five 
observation  heights  above  the  load  coil  and  at  an  RF 
input  power  of  0.75  kW  (filled  circles)  and  1.0  kW 
(open  circles) . 

Figure  10.  Electron  number  densities  at  3  mm  off  the  plasma  axis, 
at  five  observation  heights  above  the  load  coil  and  at 
an  RF  input  power  of  0.75  kW  (filled  circles)  and  1.0  kW 
(open  circles) . 

Figure  11.  Electron  temperatures  (open  triangles)  and  gas-kinetic 
temperatures  (filled  triangles)  at  3  mm  off  the  plasma 
axis  at  five  observation  heights  above  the  load  coil  and 
at  an  RF  input  power  of  1.0  kW. 

Figure  12.  Electron  number  densities  measured  directly  by  Thomson 
scattering  (filled  circles)  on  the  plasma  axis  at:  five 


iki:.,‘i  v.i:  ion  heights  .ilmVf  tlii’  lii.i'l  i-nj  !  .iml  .it  i'.F  input 

i'  «>i  !)..■'>  kW  ,  and  tlii-  i’.i  i  u  u  1 .1 1  i  ■■  1  l.TK  i- !  ic  i  !o;i  number 

densities  (open  cirolos),  obtained  iron:  the  Saha 
equation  usin|;  the  corresponding  measured  electron 
temperatures  shown  in  Figure  7. 

igure  13.  Electron  number  densities  measured  directly  by  Thomson 
scattering  (filled  circles)  at  3  mm  off  the  plasma  axis 
at  five  observation  heights  above  the  load  coil  and  at  an 
RF  input  power  of  1.0  kW,  and  the  calculated  LTE 
electron  number  densities  (open  circles),  obtained  from 
the  Saha  equation  using  the  corresponding  measured 
electron  temperatures  shown  in  Figure  11. 
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